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ABSTRACT 

Context. Coronal waves are an important aspect of the dynamics of the plasma in the corona. Wavefront dislocations are topological features 
of most waves in nature and also of magnetohydrodynamic waves. Are there dislocations in coronal waves? 

Aims. The finding and explanation of dislocations may shed light on the nature and characteristics of the propagating waves, their interaction 
in the corona and in general on the plasma dynamics. 

Methods. We positively identify dislocations in coronal waves observed by the Coronal Multi-channel Polarimeter (CoMP) as singularities in 
the Doppler shifts of emission coronal lines. We study the possible singularities that can be expected in coronal waves and try to reproduce the 
observed dislocations in terms of localization and frequency of appearance. 

Results. The observed dislocations can only be explained by the interference of a kink and a sausage wave modes propagating with different 
frequencies along the coronal magnetic field. In the plane transverse to the propagation, the cross-section of the oscillating plasma must 
be smaller than the spatial resolution, and the two waves result in net longitudinal and transverse velocity components that are mixed through 
projection onto the line of sight. Alfven waves can be responsible of the kink mode, but a magnetoacoustic sausage mode is necessary in all 
cases. Higher (flute) modes are excluded. The kink mode has a pressure amplitude that is smaller than the pressure amplitude of the 
sausage mode, though its observed velocity is larger. This concentrates dislocations on the top of the loop. 

Conclusions. To explain dislocations, any model of coronal waves must include the simultaneous propagation and interference of kink 
and sausage wave modes of comparable but different frequencies, with a sausage wave amplitude much smaller than the kink one. 
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1. Introduction 

Magnetohydrodynamic transverse waves seem to be a 
relevant constituent in the dynamics of magnetic and 
plasma structures in the solar atmosphere. Their pres¬ 
ence has been invoked to explain imaging and spectro¬ 
scopic signatures of periodic plasma motions detected in 
different types of structures, wit h different physic a l con- 


wave heating processes (iParnell & De Moorte]Ll2012tlArregui . 
201.51) . 


Lonez Ariste et al. ( 2013h demonstrated the existence of 


solutions to the equation of magneto-hydrodynamic waves car¬ 
rying wavefront di slocations. In that work, the equation studied 
was Eq. 4.14 from lPrie^ (Il982h which reads: 


Nakariakov et al 


(iDe Pontieu et al 
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ditions, such as coron al loops (lAschwanden et al.[ H999t + [V x (V x (vi x Bo))] x 


Bn 


i^Po 
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20071) : soft X-ray c oronal jets (ICirtain et al. 


2007h: pr omine nce fine structures lOkamoto et al.l (12007 ): 


(iTomczvk et al 


Lin et al.l (120091); or extended regions of the solar corona 


20071) . In recent years, their relevance has 
se of their potential as a seismology di- 


agnostic tool ( 

Arregui et al. 


.007; Goossens et al., 
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2012) and their possible role in 
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where Vi is the velocity of the plasma, Bq the magnetic field, 
assumed constant, cj the speed of sound in the medium, and 
po the density. This equation describes linear magnetohydro- 
dynamical waves in a homogeneous, isothermal medium. This 
equation has as solutions waves propagating a long the mag¬ 
netic field that carry wavefront dislocations ( Nve & Berrvi 
1974h . that is, singularities in the phase of the wave. Examples 
of such waves can be seen in Eig. [T] The four images show dif¬ 
ferent kinds of dislocations made by varying the parameters fi 
and (5 of a generic solution to the longitudinal (that is, along 
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the magnetic field ) component of the veloci ty perturbation 
inEq. ([T]l given bv Lonez Ariste et alJ ( 2013h : 


V, = A + jSe'^kiz - ct)) . (2) 

Starting from the left, the illustrated dislocations are a pure vor¬ 
tex ifi - 0), two pure edge dislocations (j3 ^ 0,6 = |) and a 
mixture of vortex and edge that can be described as a gliding 
dislocation (fi + 0,6 - 0). The axes of the plot are time and 
distance, though any other sensitive choice can be made that 
captures the topology of the wavefront under scrutiny. When 
looking at this kind of pictures it is important to agree in that 
the quasi-periodic variation seen in those time-distance plots 
is a sound representation of the phase of the wave times an 
amplitude that does not have to be constant over the plot (ac¬ 
tually it will be zero at the singularity and only at that point, 
but can take any other value elsewhere). A wave period can 
be measured with more or less precision and significance by 
measuring the distance between two crests (white color) or two 
valleys (black color) in such a plot. If we were staring at a plane 
wave, the pattern of black and white crests and valleys of the 
wave would be roughly parallel over the plot. We see that this 
is not the case and we have encircled a particular region in one 
of the figures in which one crest ends suddenly where two val¬ 
leys merge into just one. It is evident that the ending point of 
that crest cannot have a well-defined phase, it is singular. The 
amplitude of the wave at this point should be zero. This is the 
dislocation. Sometimes it is just a point (pure edge dislocation) 
but often it is a line of singularities (as in the vortex case on the 
leftmost figure or the gliding dislocation of the rightmost one). 
At this point it is important to distinguish a dislocation from 
mere wave nodes. A node is just a zero of a wave with real am¬ 
plitude. It is often illustrated in one-dimensional plots of waves, 
often two interfering waves as in standing waves, where we see 
the amplitude go to zero. But this definition makes no mention 
of the phase of the wave which, despite the zero in amplitude, 
may still be perfectly well defined. In a dislocation the phase of 
the wave is singular. Thus, dislocations can be seen as nodes, 
but not all nodes are dislocations. An illustration is offered in 
our context of MHD waves by a solution of a wave propagating 
along the z direction and which in the transverse plane (r, 0) has 
the form 7i(r)e'®. Each zero of the Bessel function J\ is a node 
of the wave; the amplitude will be zero at those places at all 
times. But the node at r = 0 is a dislocation since the phase 9 
is undefined or singular at that point. On the contrary the first 
zero at r = 3.83 is a node for which the phase is perfectly 
well defined and hence it is not a dislocation. Telling apart a 
dislocation from a node with a wave given as a real function 
of just one variable is not possible. If the wave is written as 
a complex function pe‘^ - a + ib the dislocation is found as 
that place where a - b - 0 what immediately leaves the phase 
X - arctan ^ singular and the amplitude p - -i- b^ zero. If 

the wave is seen as a function of two variables (time and one 
spatial dimension) as in our examples of Eig. [T] or Eig. |2l the 
qualitative description made above also successfully identifies 
a dislocation and separates it from other nodes. But this qualita¬ 
tive picture of a dislocation in a time-distance plot of the wave- 
front or the simple rule of finding the places where both the 
imaginary and real parts of the waves are simultaneously zero 


must be complemented with a more quantitative and mathe¬ 
matically rigorous definition of the singularity. This is done 
by drawing a closed curve, called a monodromy, along which 
we integrate an appropriate parameter. In our case the mon¬ 
odromy is computed over the phase: if we describe the wave 
as a map of pe'^, with p and x real numbers that define the 
amplitude and phase respectively at each and every point, the 
monodromy of interest is 



where C represents that closed curve that can be in our example 
the circle entouring the singularity. The monodromy is strictly 
zero when there is no singularity in the area enclosed by C, and 
it is Inm if there is a dislocation, with m the topological charge 

/~vf cir» rrn 1 or*! 7 Tl^ic /^Vio-rrT<a wt ic t/~v xT/itVt 



Fig. 1. Eour examples of dislocations in the propagation of a 
wave with the time in abscissas and a distance in ordinates. 
Erom left to right, a vortex, two edges and a gliding edge, all 
with charge 1. Around the second edge dislocation a closed 
curve has been drawn with dots and squares. The integral of 
the phase along this curve, the monodromy, is non-zero. 


Lonez Ariste et al. ( 2013h pointed to observ ations of mag- 
netoac oustic waves in the sunspot umbra by Centeno et alJ 
(l2006h where dislocations as the ones illustrated could be easily 
identified visually but also by computation of the monodromy. 
Dislocations are not extraordinary solutions of the wave equa¬ 
tion but actually quite common occurrences, as those observa¬ 
tions demonstrated. In the scenario described by Eq. ([1]) there is 
a clear axial symmetry given by the constant magnetic field. If, 
for the time being, we restrict ourselve s to waves propagating 
along this z direction it is well-known ( Wen^ 


7es prop; 

iEM 


Spruit, 


1982 ; Edwin & Roberts , 1983 ; Roberts , 1981 ) that solutions to 


this equation are given in terms of families of Bessel functions 
times an azimuthal dependence e"”®. This angle 6 correspond¬ 
ing to the cylindrical azimuth coordinate is obviously singular 
at the origin of the coordinate system r = 0. Those classic so¬ 
lutions to the Eq. O carry therefore a dislocation at r = 0 for 
all the cases with m + 0. We find here a third occurrence 
of this time to refer to the azimuthal wave number of 
the solutions to the wave equation. It is sufficient to try the 
monodromy integral over the gradient of the phase to real¬ 
ize that m is exactly the charge of the dislocation introduced 
above. And thus, following the usual naming convention, a 
wave with charge m-0 will be called a sausage mode, while 
a wave with charge m - 1 will be called a kink mode. Higher 
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values of the charge m are referred to as flute modes. Finally, 
it is interesting to notice and stress at this point that any attempt 
to describe such singularities with a finite combination of colin- 
ear plane waves or Fourier components will fail. Dislocations 
require either the full infinite series of the Fourier decompo¬ 
sition or choosing a family of solutions that already carries 
a dislocation in each or most of its components. The Bessel 
functions times the azimuthal dependence e"”® are one of such 
families and an example that we should expect to find in the 
description of magnetohydrodynamic waves in the solar atmo¬ 
sphere. It is therefore expected that dislocations are observed in 
waves in the solar atmosphere, though the common description 
of waves in terms of plane waves have resulted in overlooking 
them because plane waves cannot describe dislocations. 

In the present paper we should turn into another observation 
of magnetohydrodynamic waves, this time in the solar corona. 
The observations can be seen in Fig. |2] and they will be de¬ 
scribed in terms of dislocations in the next section. Section 3 
is a first attempt to describe the observed dislocations in terms 
of generic wave solutions that require the presence of several 
interferring waves with different frequencies and propagation 
velocities. The detailed inspection of mathematical solutions in 
Section 4 will show that although most of the waves propagat¬ 
ing in coronal tubes carry dislocations, they are not visible in 
observations like the one in Fig.|2]unless the observed Doppler 
velocity is projection of both the transverse and the longitudi¬ 
nal velocity. This last one has therefore to be non-zero, what 
implies the necessary presence of at least one magnetoacoustic 
mode in the observed waves. 


2. Dislocations observed in waves propagating 
along coronal tubes 

Figure |2] shows a plot of measured Doppler velocities as a 
function of time for a particular trajectory in the corona, 
measured in a coronal emis sion line by the Coro nal Multi¬ 


channel Polarimeter (CoMP, |Tomcz^k_et^ 


servations wer e presented by iThrelfall et al 


Tomczvk et al.L 120071 iTomczvk & McIntosh! 


2008). The ob- 


(2013) (see also 


20091) . The pa¬ 


rameters of the observation, off the solar disk, imply that we 
are observing here the velocity of the emitting plasma roughly 
transverse to the direction of the coronal magnetic field. The 
periodicity of the signal is evident. An average frequency w 
can be estimated and we should use it throughout this paper to 
describe a carrier wave whose amplitude and phase may be lo¬ 
cally modified. Another clear feature of the observed waves is 
the tilt of some of the crests and valleys, indication of a prop- 
agati ng wave along the coro nal feature and along the magnetic 
field (IThrelfall et al.L l2013l) . Several dislocations are also vis¬ 
ible, more easily after comparing them with the previous ex¬ 
amples. We have marked one of them around minute 55, as in 
Fig[T] with a closed curve made of 4 segments. A monodromy 
of interest would be the integration of the phase of the wave 
along this closed curve. Due to the importance, before proceed¬ 
ing, of classifying the observed wave pattern as a dislocation, 
we are going to compute the integral over the monodromy. Fig. 
12 shows the actual Doppler shift or velocity of the plasma, not 
its phase. We cannot therefore integrate directly the measured 


values over the curve. In the Appendix we show how the com¬ 
putation of the integral directly from the data can be made, but 
here we are going to take a more heuristic approach: One can 
always deformate the curve so that it is made of 4 segments 
along which we can safely identify the phase of the wave or its 
changes from the observations and integrate these inferred val¬ 
ues. The first such segment follows the blue valley to the left 
of the dislocation from roughly position 130 to position 50 in 
the ordinates. We will draw this segment as to follow the phase 
0 (let us agree on that valleys (blue) are at phase 0 and crests 
(red) at phase ji). Given the interpretation of these blue and red 
stripes as the valleys and crests of a wave, we should also agree 
in that roughly at the center of the referred valley there is a 
continuous line at phase 0 over which we draw the segment. A 
similar segment is drawn over the next blue valley to the right, 
also along the 0 phase value. The integral J dx along either one 
of these two segments is 0 since there is no change in phase. We 
next join the two segments with a straight horizontal line at or¬ 
dinates position 50 from left to right. This straight horizontal 
segment goes from a point at phase 0 to the next red crest at 
phase n and then ends in the next valley at phase 2ji = 0. The 
integral J dx - 27 t along this segment. Similarly, at ordinate 
point 130 we draw a horizontal and straight line from right to 
left joining the two vertical blue segments. This time the hor¬ 
izontal line starts at phase 0 and goes over two red crests and 
one blue valley before ending up in the final blue valley where 
we drew the segment. The integral is J dx = -2 x In, where 
the minus sign comes from going from right to left, rather than 
in the other sense. The full integral along this closed curve is 
therefore 


i 


dx - 2n — 4n — —2n, 


different than zero. Hence, it has a charge m - -I given the 
chosen orientation of the curve. A direct computation of the in¬ 
tegral confirms the result, as seen in the Appendix. From the 
theory of functions of complex variable, a closed integral over 
a function of complex variables is non-zero when it encloses a 
number of poles or singularities larger than the amount of ze¬ 
ros. We conclude that somewhere inside this closed curve there 
is a singularity that result in a nonzero, non trivial, monodromy. 

This computation of the monodromy can be repeated for 
many other similar points of the figure, in particular for the line 
of dislocations marked with a black tilted line that we identify 
as a possible gliding dislocation after comparison with Fig. [1] 
On the other hand, the main dislocation marked with the closed 
curve, and whose monodromy we computed, appears to be of 
the edge type. No pure vortex is seen. In order to explain these 
dislocations it is important to stress that the axes of the plot are 
time f and distance z along a coronal loop as projected onto the 
plane of the sky. The edge dislocations seen in the figure are 
therefore localized at particular values of t and z, but with free¬ 
dom about the transverse coordinates x, y or r, 6 in the more 
appropriate cyclindric reference system . Other assumed con- 
straints beyond z being along the loop (IThrelfall et al.L 1201 3l) 
are that the magnetic field follows the coronal loop and the 
wave propagates along this magnetic field. Since coronal loops 
appear as structures in coronal emission lines whose intensity 
variations are related to density enhancements, we further as- 
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All the analytical solutions found and given of magnetohy- 
drody namic wayes in coronal loopsjEdwin & Roberts, 190 
Priest, 1982; Roberts, 1981; Goossens et ab, 2012, 2 0091) can 
be written as 


lA = fir, 0)e‘ 


i(kz~cot) 


(3) 




0 


^ with no dependence of / on either z or tiir|p. This simplification 

can be justified in \he case of" homogeneous loops*^along the ° ^ 

2Qniagnetic field and constant in time, ^j^t for periods long 
compared with the period of the wave. TTiis assumption also 
makes possible to determine and fix a> and k. The a(jtuilp^p|iiip g j 
of the function f(r, 6) changes with the assumptions andphys- '' 


Fig. 2. Observation of coronal waves in the Doppler veloc¬ 
ity of a coronal emission line by COMP, from iThrelfall et al 


(120131) . The diagram has time in abscissas and the coordinate z, 
along the loop and parallel to the magnetic field, in ordinates. 
Among the many dislocations yisible, two haye been marked. 
The closed curye (that could represent the monodromy) en¬ 
closes an edge dislocation, while the straight segment follows 
a possible gliding dislocation. 

sume that the wave is propagating in a high-density cylinder. 
The phase speed of the propagating wave, finally, has been 
measured to be in the range 700 - 1000 km/s. Such speeds 
are much larger than the sound speed of the corona and are 
comparable to the Alfven and kink speeds characteristic of the 
propagation of Alfven waves and fast magnetoacoustic body 
waves in coronal loops respectively. 

3. Interpretation of the observed dislocations in 
coronal waves 

A dislocation is a phase singularity. At the exact point of the 
phase singularity, the amplitude of the complex wave has to 
be zero. This statement allows us to search for dislocations as 
zeros of the complex amplitude of the wave. We want to de¬ 
scribe a wave if/ propagating at a characteristic frequency at 
and wavenumber k along the z axis, parallel to the magnetic 
field and the coronal loop. One common form for the solution 
of such a wave is written as 

We notice that / may depend on the three cylindric coordinates 
and on time, in the most general case. Since the term 
can be neither zero nor singular it is obvious that all the infor¬ 
mation on the dislocations of this wave is contained in the / 
function that should be complex in general. We are therefore 
looking for those times and places (r, 6, z, t) where the complex 
amplitude 

||/(r,0,z,f)ll =0. 

A plane wave would have 


I f(r, 0) changes with the assumptions ancf phys 
ical phenomena considered, but it always retains these depen¬ 
dencies. If such a wave carries a dislocation, where the phase 
is singular, it will be found at those points (r, 0) where the am¬ 
plitude of the wave satisfies the condition 

ll/(T0)ll =0. 

This equation will be valid for all values of z and t. The dislo¬ 
cation therefore will be localized in the transverse plane (r, 0) 
but not in the plot (z, t) of Fig. |2] This is exactly the opposite 
situation to the one observed. The description of the wave as 
/(r, cannot produce a point dislocation in the plane 

(z, t) and cannot therefore describe Fig.|2] independently of the 
actual form of /(r, 0). In view of this and since the observed dis¬ 
locations are certainly localized in (z, f), we may doubt about 
the description of the wave as given in Fq. ©. Waves propa¬ 
gating along the field have to have this functional form as long 
as the assumptions made above on the loop hold, so we could 
as an alternative consider that the observed wave is propagat¬ 
ing across the field. This possibility was however quickly dis¬ 
carded after inspection of the possible analytical solutions to 
such a wave. We are not going to give here the details. 

Another possibility we have inspected is that what we 
are observing is the interference of more than one wave. 


Thi s possibility has been p u t forw ard by IThrelfall et al.l (120131) 
and Tomczvk & McIntosh ( 20091) from actual inspection and 
Fourier filtering of the observations. Although it is unclear 
whether the methods used by those authors are still valid in 
the presence of phase singularities, let us follow them and sug¬ 
gest the possibility of two waves propagating with different fre¬ 
quencies 0)1 and 0 ) 2 - Let us further ass ume that the frequency 
difference Aw = W 2 - wi is small, since iThrelfaU et alJ (l2013» 
filtered the observations in the Fourier space. This wave can 
be written as 


f(r, ^ gyikz-o^t+a) ^ 

(f(r, 6) + g{r, 0) - g(r, 0)i(Aajt - a)) 


(4) 


a complex amplitude multiplying the propagation term 
Assuming for simplicity, but without loss of gener¬ 
ality, that / and g are real amplitudes, we find that this wave 
can carry a dislocation when the two conditions 


with A real. This amplitude cannot be zero, unless A = 0. Plane fi^ + 8(r, 0) - 0 
waves carry no dislocations therefore. g{r, 6>)(Awf - a) = 0 


( 5 ) 
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are simultaneously satisfied. Notice that the two equations are 
the real and imaginary part of the complex amplitude of the 
wave. If we require that both are simultaneously zero, we get 
that the real amplitude of the wave is zero and that the phase 
is undefined, a dislocation. The second one of those equa¬ 
tions, the one for the imaginary part, is immediately satisfied 
at f = -alAco (modulo n) and so this wave may carry a dis¬ 
location localized in time at the position (r, 6) where the real 
part of the complex amplitude is also zero. It is sufficient that 
two waves, even plane waves, with slightly different frequen¬ 
cies interfere for a dislocation localized in time to be possible. 
In particular two waves with same amplitude and in antiphase 
f - -g will always carry one dislocation at time t - a/Aw. 
Such dislocation would appear in observations like the one of 
Fig.|2]as a vertical line of dislocations. We do not observe such 
dislocation. We further need to localize the dislocation at a par¬ 
ticular position z along the coronal loop. 

A first attempt would be to suppose that the two interferring 
waves also have slightly different wavenumbers k. It can be 
easily seen that one of the conditions for a dislocation in such 
a case would be 

ig{r, 0)(-Akz + Aojt - a) -0 (6) 


4. Simulating observed dislocations 

Since magnetoacoustic waves with different velocities could 


□ 


□ 



obs 


The solution to this equation is a straight line z - of dis¬ 
locations with slope the ratio of differences in wavenumber 
and frequency. This is a encouraging result, since it could be 
the explanation of the tilted gliding dislocation marked with a 
black line in Fig.|2] Although in our explanation the line of dis¬ 
locations is continuous, while in the observation is limited to 
roughly 3 periods, we could propose that we are seeing there 
an interference of two waves with slightly different wavenum¬ 
ber and frequency, the ratio of which can be measured in the 
slope of the line, resulting in a dislocation that surfs the wave 
along the coronal loop. 

But such cannot be the explanation for the more common 
edge dislocations observed. We return to the general expression 
of Eq. (O and we see that if f{r, 6) is complex, the condition 
of zero amplitude translates into both the real and imaginary 
parts being zero independently. If we require a dislocation to 
be localized in two coordinates, like z and t, both the imagi¬ 
nary and real parts have to depend on them. The interference of 
two waves with similar wavenumber and/or frequency allows 
us to introduce an imaginary term ig(r, 0){-Akz + Aa>t) to the 
complex amplitude of the wave, but not to the real part. We end 
up therefore with either just one coordinate fixed for the dis¬ 
location or with a linear relationship between both, as we saw, 
but not with a complete determination of both coordinates. To 
achieve this we will need the condition on the cancellation of 
the real part of the complex amplitude to involve z, t or both. In 
the interference of the two waves proposed in this section, the 
other vanishing condition reads 

fir, 0) + gir, 0) = 0, 

which does not depend on either z or t. The scenario of two 
interferring waves can at most explain one kind of observed 
dislocation, the gliding edge, but not the most common one, 
the edge dislocation, in the observations of Fig. |2] 


Fig. 3. Cartoon showing the loop, the transverse and longitudi¬ 
nal velocities of the waves and their compination into the ob¬ 
served velocity. 


A first attempt can be made with two waves having the 
same value of m. This does not work. No dislocation can be 
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observed in this configuration. The second attempt concerns 
two wave s of different value of m . This was already sug¬ 
gested by iLonez Ariste et al.l (120 13h as a means to produce 


observed dislocations. Let us suppose first that both waves 
are magnetoacoustic ones. We can suppose here a mode with 
m - 1, a kink, superposed with a mode with m = 0, a 
sausage. This comes ha ndy because the dispersion relations 
((Edwin & Robertslll983h show that these two modes propagate 
at different speeds and frequencies. The simultaneous existence 
of a sausage and a kink mode in the presence of a coronal den¬ 
sity tube will imply that, either the tube has a large enough ra¬ 
dius, or that the sausage mode is a slow mode, two alternatives 
to which we should return later. To proceed with our descrip¬ 
tion we need to have explicit expressions of the velocities of 
those wave modes. To write them we are going to place our- 
selved in a simple model with a coronal loop formed hy a 
uniform cylindric tube of plasma with piece-wise constant 
density. In the interior of that cylinder, the magnetoacoustic 
kink mode can he written as 


Vr = -Al 


ojf - k\ci 


2 2 
oj\m^ 


Jo(mir) - 

niir 


Vg - -iAi- 




a>pn^ r 




Vr - -Aq- 
Vg = 0, 


to: 


2 - 

2 2 
ajvnQ 




The transverse velocity in both cases has yet to be projected 
onto the plane that contains the line of sight and the z direction. 
For this we have first to define the direction 6-0. Without loss 
of generality, we set this direction along the line of sight. This 
choice simplifies the expressions above and for the magnetoa¬ 
coustic kink we can write 


of reducing the number of long intermediate expressions, we 
are going to introduce the last ingredient of our model. The 
interference of two waves with different frequencies fixes the 
time of occurrence of the dislocation. The position z of the dis¬ 
location is given by the position along the loop at which the lon¬ 
gitudinal and the transverse components of both wave modes 
are projected with the right angle. But, with the expressions we 
have at hand at this point, we realize that diferent points in the 
plane (r, 6) place the dislocation at different values of z and t. 
To find a dislocation at z and t independently of r and 6, we 
are going to assume that the cross-section of the coronal loop 
is smaller than the spatial resolution of the CoMP instrument. 
Coronal loops are mostly unresolved with present instruments, 
and certainly they are so with CoMP. We are therefore assum¬ 
ing that the full transverse plane (r, 6) of the wave is contained 
in one pixel, that is, that the radius R of the coronal tube is 
smaller than the pixel, and this forces us to integrate the veloc¬ 
ities in both variables r and 6: 


Vobs 


1 ^2?! 
—T I dr I I 
rtR^ Jo Jo 


rdOvio 


(7) 


where nii is the mo defined by Edwin & RobertsI ( 1983h for the 
case of the kink mode frequency and wavenumber. Always us¬ 
ing the same loop model, and with a similar redefinition of 
mo, the magnetoacoustic sausage (m = 0) mode can be written 
as 

r^kn 

Vj = -iAo-^Jo(inor)e‘’^^~“^'^‘ 


The integral in 0 is particularly interesting since it cancels out 
most of the terms in the expressions of the velocity, symmetric 
in azimuth. This cancellation of unresolved velocities has also 
been pointed as a concern for the observation of pure Alfven 
waves with m - 0, which are torsional azimuthally symmet¬ 
ric waves. While the sausage mode has to be magnetoacous¬ 
tic, we notice at this point that the kink mode could either be 
magnetoacoustic or Alfven, both cases resulting in a non-zero 
transverse velocity after integration on r and 6. In what follows 
we will pursue the calculations for the magnetoacoustic kink 
mode, keeping in mind that the magnetoacoustic sausage mode 
could also be interferring with an Alfven kink mode. 

After integration on r and 6, the observed velocity from the 
magnetoacoustic kink mode, with both the longitudinal and the 
transverse velocities combined, is 


( 8 ) 


1 or - k^.c] . 

VobsO = ^Ai-.^os// 

R to.m\ 


while the sausage mode is seen as 
2 .. koc^ 


f 


rJo(m\r)dr 


(9) 


VoKO ='"“'sin/i f rJo(mor)dr (10) 

« Wa Jo 


Vtrans — Al - 


(X> 


7 - yy 


2 2 


Finally, to combine the two waves we rewrite to - u)\ and 
0)0 - coi + Aco. For simplicity, we will assume that ko - k\ - 
'"i' y-m\Jo{m\r){cos^ 9 + isin0cos0)-i-^ since this does not alter the results. The observed velocity, 

following our model, will be 


H—7i(mir)(cos20 + i sin 20) 
r 

while in the case of the magnetoacoustic sausage the transverse 
velocity is just 


Vobs = 


1 


^trans — -^ 0 " 


coi - kic 


•2^2 


(olml 


“^o'Jdmor)cos0 


2i 


+ 2 


R-^ 

kcl 


-2^2 


k^C 


to^mi 


'f 


AiCOSyU I rJo{m\r)dr— 


r 

An sin u cos Awf | rJo(mor)dr+ 
{to + Aw)2 Jo 

kc^ r* 

--y^Ao sin u sin Awf I rJo(mor)dr 

(w -H AojY Jo 


( 11 ) 


Before combining the transverse and longitudinal components 
projected onto the line of sight for both modes, and in the sake 


This observed velocity wave will show a dislocation when both 
the imaginary and real parts of the amplitude are simultane- 
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ously zero. This leads to the following two equations 


2^2 


- k^C 
aP-mi 


AiCOSyU I rJo(mir)dr + 


kcz 


+2-^— -Aq sin 11 sin Awf I rJn(mQr)dr - 0 

(U) + A(jl)Y ' 


■I 

'i: 


kcj 


(to + AtoY 


sin yu cos Aojt - 0 


( 12 ) 

(13) 


The second of these two conditions implies that, other than the 
trivial case ji -0, the observed velocity will present a disloca¬ 
tion only if 

cos Atot - 0 

that is, at those times when the two waves are in anti-phas^B 
This fixes, as expected, a time t when the dislocation is possible 
depending on the difference in frequency and on the precise 
moment when either wave was excited at, let us say, the feet of 
the loop. It is interesting to notice that, given the wavelength 
of the coronal waves, one does not expect this cophasing to 
happen more than once per loop. By inserting this cophasing 
condition in the first condition for the dislocation we obtain 
a formula for the angle yU at which such dislocation would be 
visible; 

lAyluY 1 Aoj'^ rjQ(mir)dr 

tanju =- -\^-k^\ - 1-H- ^ -,(14) 

2 Ao \ Cj j m\k\ to I rJo(m(jr)dr 


with an equivalent formula for the case of an Alfven kink wave 
interferring with the magnetoacoustic sausage mode. 

A dislocation will be seen at time f at the position z where 
the two waves happen to be in phase and which is seen by the 
observer at an angle ji given by Eq. (fT4l) . Fig.|4]shows examples 
of the variation of the angle with the ratio of amplitudes for 3 
different values of the phase speed of the magnetoacoustic kink 
wave. It has been computed for typical values of the sound and 
Alfven speed in the corona (100 and lOOOkm/s respectively) for 
an average wave period of 3 minutes and a 5% difference in fre¬ 
quency between the two waves. As an example of what would 
be seen, we show in Fig. |5]the observed Doppler velocity for 
the case of a fast kink wave propagating at 7 times the speed of 
sound (700 km/s), interferring with a slow sausage wave, that 
presents a ratio of amplitudes ^ = 0.25. Notice that these are 
the scalar amplitudes of the full vector wave or, in other words, 
the amplitudes of the pressure wave. The observed velocities 
have amplitudes given by these Ai and Aq times other factors 
(see Eqs.|7]and[8]l resulting in the kink velocity having a larger 
amplitude than the sausage velocity even if Ai is smaller than 
Aq. With these numbers, as expected, a dislocation appears at 
roughly ji - 50°. The result is visually striking as a correct re¬ 
production of the edge dislocations observed in coronal waves 
by CoMP and seen in Fig. 

If sausage and kink waves were excited in perfect phase 
matching at the loop feet, the dislocation would be visible peri¬ 
odically at most at one position along the loop (with a very long 
period P - 2^). The observed changes in position and time 

* Had we used the Alfven kink (m = 1) mode instead, the condition 
here would have been that both waves have to be in phase. 


mean that the conditions of excitation of kinks and sausages in 
the loop feet vary permanently. Notice that, if the kink wave 
propagated at roughly the speed of sound, the factor ^ - k\ 

\x/rvnlH cinH Hiclrvr*i:itir»n \x/rviilH rvnl\7 f/AiinH iit lrv\x/ 


□ 



Fig. 4. Angle yu at which the dislocation will be seen following 
Eq. (fT4l) for different velocities of the magnetoacoustic kink 
wave, assuming a sound speed Cj = IQQkmfs, an Alfven speed 
10 times the speed of sound, a period of 3 minutes for the kink 
wave and a 5% difference in frequency between the kink and 
the sausage waves. 


5. Discussion 

Observations of waves in the Doppler measurements of coronal 
emission lines by CoMP show wavefront dislocations. Several 
qualitative and quantitative arguments have been made to en¬ 
sure this identification of the wave singularity in the observed 
data. The conclusive identification of these singularities called 
dislocations requires an explanation in terms of the waves ex- 

















3 . 



Fig. 5. Simulation of an observation of magnetoacoustic kink 
and sausage modes propagating along a circular loop with a 
5% difference in frequency, at a ratio of 7 in phase speeds and 
a ratio of amplitudes ^ = 0.25. The z coordinate along the 
loop is given in terms of the angle yU it forms with the line of 
sight at each point. 


pected to propagate in those coronal regions and their velocities 
projected along the line of sight, seen as Doppler shifts in the 
coronal emission lines observed by CoMP. 

We have provided such an explanation while explicitly dis¬ 
carding several other possibilities. We recall that the observa¬ 
tions were projected in a diagram z—t where z was the direction 
of propagation of the waves, also assumed to be the direction 
of the magnetic field. The dislocations were mostly of the edge 
type, meaning that the phase singularity was found at isolated 
points in the z—t plane. The waves naturally propagating along 
the magnetic field in the corona carry dislocations, but those 
dislocations are of the edge type only on the plane transverse 
to z, which in cylindric coordinates we can refer to as the plane 
r-0. These basic dislocations carried by the waves propagating 
along the magnetic held in coronal conditions cannot explain 
the observations. Our hrst conclusion in this paper is therefore 
that the observations cannot be explained by single propagat¬ 
ing waves of one type or the other: the presence of dislocations 
forces us to look for more elaborated scenarios of propagating 
waves. 

Three ingredients are necessary to hx the dislocation in one 
single point of the z-t diagram as observed. First, the interfer¬ 
ence of two waves with different frequencies can hx the dislo¬ 
cation in time, but not in z. Second, the combination of trans¬ 
verse and longitudinal velocity components by projection onto 
the line of sight, can hx the dislocation in z but not in time, as 
long as the transverse and longitudinal components belong to 
waves of different mode m (sausage and kink waves for exam¬ 
ple). Furthermore, the position z of the dislocation varies for 
different points in the transverse plane r-6. Third, the integra¬ 
tion of the line emission over the transverse plane r-9 allows to 
hx a single z and t for the whole emitting region with two inter- 



ferring waves at different frequencies, if that transvefce sHtio 
of the loop is not spatially resolved. 

In our analysis, we used quite a simple model | 
nal loop and the propagation of waves]along 
phisticated or realistic models, adding |in parti| 
homogeneous density profiles across the loop r^ 
the piecewise-constant density we used, will sti 
the azimuthal symmetries in^the solutHihs to the 
ing modes that are at the core of our explaiution(of the oh 
■served dislocation^ The actual detailsj|of the s&tion, tl 
radial dependencies in terms |of l^ssel function^r other 
■he distr||hutibn offvisihle dfslocationllalong tn loop, a 
ihosVaspecis mayLchange in ihosCf^ore realisne models! 

It is ■Disturbing thni; for exampldf when the derisiiy is al- 
^idn^f^TO^^^ntinuous manner, the classic kink wave 
beco mes a surface Alfven wave w ith a strong Alfvenic char¬ 
acter [CoossenTet^ (I2009il20l3) and the solutions display 
rapid variations in the radial and azimuthal components 
at the resonant layer with the consequence that the global 
motion is quickly damped. It is difficult to foresee what the 
consequences would be for the existence and properties of 
dislocations as described in this study. But we are confi¬ 
dent that the need for interference of two modes with dif¬ 
ferent azimuthal number m at slightly different frequencies 
with longitudinal and transverse velocities combined in the 
projection onto the fine of sight will still be the ingredients 
to explain the observed dislocations. Simpler conditions do 
not seem to work, and more complicated ones may not be 
so common and harder to be simultaneously met. 

Indeed, we appreciate in our explanation that the three 
ingredients (interference of waves of different frequency, com¬ 
bination of longitudinal and transverse velocity components 
into the line of sight, and integration of the full emitting re¬ 
gion inside the pixel) are quite natural and common. It is not 
surprising therefore that the observations show a large number 
of dislocations during the observation. But the same three in¬ 
gredients also limit the type of waves responsible for the obser¬ 
vations. Thus the integration over the r-9 plane excludes from 
the model any wave with azimuthal symmetry. Pure Alfven 
waves with m = 0, in particular, result in a zero signal and 
can be excluded. The need for a longitudinal velocity compo¬ 
nent implies the presence of at least one magnetoacoustic wave. 
Sausage modes (m - 0) can be responsible of the longitudinal 
velocity component, but their transverse velocity is azimuthally 
symmetric and cancels out. Kink waves (m — 1), either mag¬ 
netoacoustic or Alfven, can be responsible of the transverse 
velocity component, but not of the longitudinal one, in one 
case because it cancels out, and in the other because there is 
no longitudinal component. The addition of a magnetoacoustic 
sausage plus a kink wave, either magnetoacoustic or Alfven, 
with slightly different frequencies, has all the required proper¬ 
ties. Integrated across the loop, the observed Doppler shift is 
made of the transverse velocity of the kink mode plus the lon¬ 
gitudinal velocity of the sausage mode. At those points along 
the loop with the right projection angle, the interference of the 
two waves with the right phases results in a dislocation local¬ 
ized in z - f as observed. The equation for the appearance of 
the dislocation also suggests that, at high speeds of the kink 
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wave, dislocations will be more frequently observed on the top 
of the loops, but some variability in position is expected if the 
amplitude of the kink mode is smaller than that of the sausage 
mode. All this appears to be the case in the solar corona, and 
therefore we conclude that a fast kink mode, either a fast mag¬ 
netoacoustic or and Alfvn one, catches up a slow sausage mode 
at some point along the loop and produces a dislocation visible 
with COMP. The fast character of the kink mo de is suggested 
by th e measured speed of the observed wave ( Threlfall et al.L 


20131) . but also by the comparison of the predictions of Fig. |4] 


with the observed positions of the dislocations in Fig. |2] The 
slow character of the sausage mode makes it compatible with a 
propagation with cross-sections smaller than CoMP spatial res¬ 
olutions that justify our integration in r and 6, and apparently 
smaller than the cutoff of fast sausage modes. A slow sausage 
mode makes it also a good candidate to be over-run often by 
fas t and/or Alfven kink m odes excited at different times. 


Threlfall et al.l (12013h compared the observed waves in 


CoMP with simultaneous observations of emission, and hence 
density, perturbations observed by AIA. The density perturba¬ 
tion corresponding to our two waves is 


0) 


(15) 


Despite the absence of any azimuthal dependence, the radial 
integral of the Ji(mor) in the case of a fast magnetoacoustic 
kink is going to be almost negligible compared to the same in¬ 
tegral of the Joimor) function for the slow sausage mode. The 
compressional wave will therefore be dominated by the sausage 
mode. The Doppler signal, on the other hand, will be dominated 
by the kink mode. Hence, the two instruments are sensitive to 
one or the other wave but not to both, and one should not expect 
any spatial correlation between observations of the compres¬ 
sion wave by AIA/SDO and observations of the velocity wave 
by COMP. In spite of this, the observed periods (frequencies) 
w ill of course be simila r. This coincides with the conclusions 
of lThrelfalleTaD(l2013h . 


6. Conclusion 

We have identihed wavefront dislocations in the the observa¬ 
tions of coronal waves made by CoMP. Explaining the ob¬ 
served dislocations has forced us to abandon the image of a sin¬ 
gle wave propagating in those coronal structures. We explain in 
detail the only scenario we have found that can explain the ob¬ 
servations. Our model is made of two propagating waves with 
different wave frequencies and in interference. The two waves 
have different azimuth dependence (or charge) and the obser¬ 
vations integrate the velocity wave all over the cross-section of 
the wave, smaller than the spatial resolution of the instrument. 
This eliminates many possible candidates for which the sig¬ 
nals cancel out after integration. In particular torsional Alfven 
waves with m - 0 are excluded but magnetoacoustic sausage 
waves appear as a necessity, combined with a kink mode which 
can be either a fast magnetoacoustic mode or an Alfven wave. 
The two wave modes, the magnetoacoustic sausage and the 
kink, propagating at different frequencies and integrated over 
the loop cross-section are seen under different projection an¬ 
gles at different positions along the loop. This projection of the 


transverse and longitudinal velocities of the two waves onto the 
line-of-sight fixes when and where the dislocation will be seen. 

Computation of these conditions for the visibility of the dis¬ 
locations leads us to conclude that our model can reproduce the 
observations if we assume a fast kink mode (magnetoacoustic 
or Alfven) catching up over a slow sausage mode propagating 
along the coronal loop. The observed signature is dominated by 
the velocity amplitude of the fast kink mode although the pres¬ 
sure amplitude of the slow sausage is still larger than that of 
the kink mode. The observed dislocations also imply, follow¬ 
ing our model, that the spatial resolution of the observations 
is not enough to resolve the cross-section of the loop; we see 
an integrated signal. Despite the simplified model used for 
our analysis, the conditions under which dislocations can 
he seen in the data appear quite general and translatable to 
more sophisticated models. 

These sausage modes required in our scenario to explain the 
dislocation may coincide with th e density waves observed by 
AIA/SDO ( Threlfall et al. . 2013h . The density perturbation in 
our model is dominated by the sausage mode, even if the kink 
mode is magnetoacoustic. Observations of emission changes, 
mostly proportional to the density, would then mostly see this 
mode, rather than the kink one. Although of similar frequency 
and period, the very different wavelengths will make them diffi¬ 
cult to correlate with the observed velocity wave, mostly dom¬ 
inated by the kink mode. 

Finally, we insist in that the observations cannot be ex¬ 
plained with Alfven waves alone, but require a combination of 
at least one magnetoacoustic mode. The relative amplitude of 
these modes will have to be explained by the excitation of these 
waves or its propagation and eventual dissipation in the low 
parts of the corona. The existence of other propagating modes 
cannot be fully discarded, but since they cannot be responsible 
for the observations due, as for the case of the torsional Alfven 
wave, to their azimuthal symmetries, they have to be consid¬ 
ered as non-observed or without relevance for the present ob¬ 
servations. On the other hand, the continuous existence and in¬ 
teraction of two waves must be seen as a permanent feature of 
coronal loops. 


Appendix A: Direct computation of the 
monodromy 

In the main body of this paper we computed the monodromy 
through the smart trick of identifying in the data places of 
known phase and drawing the closed curve through those 
places. We drew lines along contiguous places of known phase 
in the crest or the valleys of the waves at both sides of the 
suspected dislocation and joined them with straight line that 
crossed an integer amount of crest or valleys. This method was 
suggested by the organized patterns in the data that allowed an 
easy and safe identification of those places of known phase. 
But in order to feel fully confident about the method and the 
presence of the dislocation one would prefer to be able to di¬ 
rectly solve the monodromy integral along any closed curve in 
a general wave field. 

The observed data in our present case is the Doppler shift 
of an emission line interpreted as the line-of-sight velocity of 
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the plasma. This is a real quantity. Our first step will be to inter¬ 
pret these observations as the real part of a complex wave field 
of which we have to determine the imaginary part. Let us for 
simplicity assume that the observation can be safely interpreted 
as due to a wave with a unique average frequency a>. We can 
describe the observed wave field as follows 


^(z, t) - A(z, t) cos (cut -H a(z, t)) 

The observations, the real quantity (p at each position (z, t), 
are described as a variable real amplitude A{z, t) times a co¬ 
sine variation in time with frequency w. We asume a constant 
zero time for the full wave field, but allow for a local phase 
shift a(z, t) which, through its time dependence, may include 
local frequency changes. The combination of the variable am¬ 
plitude and local phase shifts allows the description of very 
complicated wave patterns, including the one in Fig.|2] as long 
as one accepts the constant average frequency co over the time 
and place of the observation. 

We can decompose the cosine function as 

cos (wf + a(z, 0) = cos a>t cos a - sin ojt sin a 


The local phase shift a can now be interpreted as a local modi¬ 
fication of the amplitude of two different waves: 

0 = [A{z, t) cos a cos wf] - [A(z, t) sin a sin ut\ - 

- tpc cos ojt + ips sin tot 

This suggests the construction of the complex wave field 
piz, f) = {pc + 

This wave has an amplitude 

+ A^ 

the amplitude of the observed wave, while locally it adds a 
phase 

ips sin a 

tan y = — = -= tan a 

ipc cos a 

identical to the local phase of the observed wave. Thus the pro¬ 
posed complex wave field 


has the same observable parameters as the original real wave 
4>{z, t) and can be used instead of it, with a straightforward (dif- 
feomorphic) correspondence between them. 

Since z and t are the coordinates of Fig. |2] let us re-write 
this complex field for the wave on the longitudinal velocity as 
just a real amplitude and a phase 

ip(z, t) = pe'^. 


After differentiating this expression, we find that 
dip - dp e'^ + ip e'^ dx- 


Dividing by ip 

dip dp 
— = — + idx . 

p p 


(A.l) 


(A.2) 


We can integrate both sides of this last expression along a 
closed curve C, the monodromy: 



dx- 


Since p, the amplitude of the complex wave, is, by definition, a 
real quantity we have that ^ y = 0 on any closed curve, and 
the monodromy simplifies to 



(A.3) 


The right part is the monodromy over the phase of the wave 
which, if different than zero, identifies the presence of a sin¬ 
gularity, a dislocation, inside the closed curve. The left part is 
an integral over the observed data. We conclude that from the 
observations we can built de complex field p and then compute 
the integral on the left along the chosen closed path C to obtain 
the required integral over the phase at right. This solves the 
problem of computing the monodromy on the phase directly 
from the data. It is illustrating to make one further step. The 
integral on the left side can be formally integrated 


'jPl 



(A.4) 


where pf and pi are the final and initial values respectively of 
p at the closed path. These would be the same, since the path 
is closed, but for the p function being complex. This opens the 
possibility of those initial and final values not being in the same 
Riemann surface of the logarithm function. Indeed the complex 
logarithir0 is log z = (log zjpUncipai + i2nN. The principal part of 
the logarithm is identical for the initial and final points, so that 
we can conclude that 



(A.5) 


where N is the number of twists made by the logarithm function 
as it follows the path C. 
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